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1.  Introduction 


Antimony  triiodide  (SbE)  is  one  of  the  most  effective  n-type  dopants  for 
(Bi,Sb)2(Te,Se)3  thermoelectric  materials.1'2  Precise  control  of  the  concentration  of 
the  Sbp  dopant  in  (Bi,Sb)2(Te,Se)3  is  the  key  to  obtaining  good  n-type  materials, 
and  consequently,  the  highest-performance  thermoelectric  devices.  There  is  an 
ideal  amount  of  SbE  that  yields  an  optimum  concentration  range  of  3  to  5xlOly/cm3 
of  n-type  carriers.2  All  commercial  bulk  thermoelectric  devices  have  doped  n-type 
materials  prepared  to  achieve  this  n-type  dopant  concentration  range. 

To  effect  n-type  doping  in  bulk  materials,  an  appropriate  amount  of  SbE  is  added 
to  (Bi,Sb)2(Te,Se)3  source  materials,  which  are  all  co-melted  at  temperatures 
greater  than  600  °C.  In  the  melt,  SbE  chemically  decomposes  and  liberates  free 
antimony  (Sb)  and  iodine  (I)  atoms,  which  become  substitutionally  incorporated 
onto  the  (Bi,Sb)2(Te,Se)3  crystal  lattice  during  subsequent  solidification  and 
Bridgman  crystal  growth.  Substitutional  Sb  is  isoelectronic  on  a  bismuth  (Bi)  site 
and  has  little  consequence  on  doping.  However,  substitutional  I  on  a  tellurium  (Te) 
site  will  serve  as  a  singly  ionized  donor,  and  in  the  case  of  I  on  a  Bi  or  Sb  site,  will 
act  as  a  doubly  ionized  donor.  In  either  of  those  cases,  substitutional  I  causes  n-type 
conduction. 

The  growth  of  epitaxial  single-crystal  thin  films  of  n-type  (Bi,Sb)2(Te,Se)3 
materials  presents  new  doping  challenges  because  it  is  a  nonequilibrium  process. 
(Bi,Sb)2(Te,Se)3  epitaxial  materials  have  been  grown  by  molecular  beam  epitaxy 
(MBE)3  and  organometallic  vapor  phase  epitaxy  (OMVPE)  growth  techniques.4 
However,  the  use  of  this  preferred  SbE  dopant  in  either  MBE  or  OMVPE  is  not 
reported.  MBE  and  OMVPE  are  distinct  from  the  bulk  technique  because  crystal 
growth  proceeds  by  the  vapor  phase  and  at  a  significantly  lower  temperature, 
approximately  300  °C.  Because  the  material  never  melts,  the  incorporation  of  SbE 
into  epitaxial  films  would  necessarily  proceed  by  some  alternate  vapor  phase 
mechanism.  The  goal  of  the  present  work,  therefore,  is  to  understand  the  vapor 
pressure  behavior  of  SbE  as  a  stepping-stone  toward  the  overall  goal  of  employing 
SbE  during  MBE  growth  to  obtain  extrinsically  doped  n-type  (Bi,Sb)2(Te,Se)3 
epitaxial  materials. 

2.  Vapor  Pressure 

In  a  classical  MBE  growth  process,  each  individual  source  material  is  heated  to 
some  temperature  (T)  to  produce  an  equilibrium  vapor  having  pressure  (P)  that  can 
be  made  to  condense  on  a  crystalline  substrate  and  epitaxially  grow  single  crystals 
of  the  desired  chemistry.  The  production  of  an  equilibrium  vapor  pressure  above 
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each  condensed  material  is  a  phase  transformation  process  where  there  is  a  net 
change  in  molar- specific  volume  (AV)  from  the  condensed  phase  to  its  vapor.  The 
rate  of  that  phase  change  depends  on  the  enthalpy  of  vaporization  (AHyapor)  and  is 
governed  by  thermodynamics.  The  Clausius-Clapeyron  relationship  shown  in 
Eq.  1  describes  the  thermodynamics  that  governs  the  phase  transformation  between 
phases: 


dP  _  (AHvap0r)  l 

dT  ~~  AV  Lt.  ' 


(1) 


For  a  phase  transformation  from  a  condensed  phase  to  a  vapor  phase,  one  may 
assume  the  relative  change  in  molar- specific  volume  (AV)  is  essentially  equal  to 
the  specific  volume  of  the  vapor  (VVaPor)  because  vapor  phases  usually  occupy 
significantly  larger  volumes  than  those  occupied  by  condensed  phases  (Vcondensed). 


AV  =  (V, 


vapor 


V, 


condensed  j 


V 


vapor  • 


(2) 


If  one  assumes  the  obtained  vapor  behaves  according  to  the  ideal  gas  law,  then 
V vapor  =  RT/P,  where  R  is  the  common  gas  constant  (8.314  J/mol-K),  and  Eq.  1  may 
be  rewritten  as 


dP  P(AHvap  or) 

dT  ~~  RT2 


Rearranging  Eq.  3  and  integrating  between  2  thermodynamic  states,  say  state  “1” 
and  state  “2”,  as  per  Eq.  4, 

rP2  dP  _  (AHvap0r)  ~T2  dT 

Jpi  —  ~  r  Jt  1  72  (  ) 


yields  the  analytical,  closed-form  Clausius-Clapeyron  equation  for  equilibrium 
vapor  pressure  as  a  function  of  T  evaluated  between  states  1  and  2: 

In  (P)  =  42=2  or  p  =  exp  [42=2]  +  c.  (5) 

Equation  5  predicts  an  exponential  dependence  of  equilibrium  vapor  pressure  on 
inverse-temperature  with  the  slope  of  (AHVapor/R)-  One  method  of  experimentally 
determining  AHyapor-  therefore,  is  to  measure  the  equilibrium  vapor  pressure  of  any 
material  as  a  function  of  inverse-temperature  where  the  slope  of  the  data  can  be 
used  to  directly  determine  AHvapor- 
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3.  Experiment 


An  ampoule  of  99.999%  pure  SbE  was  loaded  into  a  Knudsen  cell  fabricated  from 
pyrolytic  boron  nitride  (p-BN).  p-BN  is  chemically  nonreactive  and  inert  over  the 
entire  temperature  range  of  concern,  so  there  are  no  chemical  effects  other  than  the 
phase  change  of  SbE  within.  Care  was  taken  to  minimize  exposure  of  SbE  to  lab  air 
because  it  is  hygroscopic  and  loses  its  brilliant  orange  color  (Fig.  1,  right)  as  it 
adsorbs  water  from  the  lab  air. 


Fig.  1  (left)  Knudsen  cell  showing  the  orifice  and  (right)  the  ampoule  of  high-purity  Sbl3 

A  Knudsen  cell  is  a  special,  almost  fully  enclosed  chamber  that  has  a  small  orifice 
on  one  end.  A  material  placed  inside  the  Knudsen  cell  is  heated  and  will  respond 
by  generating  an  equilibrium  vapor  pressure  inside  according  to  the  Clausius- 
Clapeyron  equation.  Because  some  vapor  can  escape  from  the  small  orifice,  the 
equilibrium  vapor  pressure  can  be  determined  by  techniques  that  can  quantify  the 
efflux  of  vapor.  Further,  the  enthalpy  of  vaporization,  AHyapor,  can  be  determined 
by  measuring  the  temperature  dependence  of  that  pressure. 

The  Knudsen  cell  assembly  was  mounted  onto  an  MBE  system  that  can  achieve 
vacuum  on  the  order  of  10~10  Torr.  Because  of  that  extremely  strong  vacuum,  vapor 
pressure  measurements  can  be  expanded  by  more  than  10  orders  of  magnitude 
compared  to  previous  work.  The  equilibrium  vapor  pressure  within  the  MBE  is 
measured  by  the  use  of  a  common  hot-cathode  Bayard- Alpert  gauge. 

For  the  present  study,  the  range  of  temperatures  where  the  vapor  pressure  is  small 
was  investigated  and  compared  to  previously  published  data  from  the  temperature 
range  where  the  vapor  pressure  is  relatively  high.  Three  independent  measurements 
of  the  vapor  pressure  were  conducted  in  this  study  and  are  color-coded  red,  blue, 
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and  purple  in  Fig.  2.  For  each  measurement,  the  temperature  of  the  Knudsen  cell 
was  very  slowly  increased  starting  from  14  °C  and  allowing  sufficient  time  for  the 
SbF  within  the  Knudsen  cell  to  reach  thermal  equilibrium. 


1/T  (Kelvin1) 


0.0010  0.0015  0.0020  0.0025  0.0030  0.0035  0.0040 


Fig.  2  Equilibrium  vapor  pressure  of  Sbl3  as  a  function  of  inverse  temperature 

For  the  analysis,  the  pressure  was  determined  as  the  average  of  the  3  separate  data 
collections  ranging  from  below  room  temperature,  -14  °C,  to  slightly  above  room 
temperature,  43  °C.  As  shown  in  Fig.  2,  the  obtained  data  were  self-consistent  with 
little  divergence.  In  this  temperature  range,  the  slope  of  the  averaged  pressure  data 
from  3  collections  is  (AHvapor/R)  and  follows  the  exponential  of  (-15,452.91/T) 
and  yields  AHvapor  =  -128.48  kJ/mol. 

Further,  the  extrapolation  to  higher-temperature  data  show  reasonable  consistency 
with  previously  published  data;  the  data  would  connect  to  form  a  piecewise 
continuous  curve.  There  is  some  nonlinearity  and  softening  in  the  global  slope  as 
the  melting  point  is  approached,  and  then  surpassed.  On  Fig.  2,  the  previously 
published  data  suggest  a  softer  (AHvapor/R)  slope,  which  is  found  to  be 
(-9719.5/T),  and  that  suggests  a  softening  of  the  enthalpy  of  vaporization  to 


Approved  for  public  release;  distribution  is  unlimited. 

4 


AHyapor  =  80.81  kJ/mol.  In  the  mid-temperature  range,  as  the  melting  point  is 
approached,  the  (AHVapor/R)  slope  moderates  to  an  intermediate  value  of 
AHyapor  =  87.6  kJ/mol.  This  progressive  reduction  is  most  likely  associated  with 
increased  translational,  rotational,  and  vibrational  entropies  and  larger- amplitude 
intermolecular  oscillation  as  the  solid  approaches  the  melting  point. 

The  melting  point  of  SbH  is  170.6  °C  and  is  indicated  on  Fig.  2  by  the  vertical  black 
line  intersecting  the  x-axis  at  approximately  1/T  =  0.00225/K.  Above  the  melting 
point,  there  are  more  divergent  data  sets  from  previously  published  work  that  imply 
the  thermal  decomposition  of  SbB  into  unknown  subspecies,  which  are  not  relevant 
to  the  present  work. 

4.  Discussion  and  Measurements 

Bruner  and  Corbett5  published  results  from  their  study  of  the  vapor  pressure  of  SbF 
at  very  high  temperatures  as  a  step  toward  understanding  the  high-temperature 
behavior  of  new  linear-molecule  structural  forms  of  halide  given  by  the  reaction 
4Sbl3  +  2Sb  =  3Sb2l4.  For  pure  liquid  Sb I3,  however,  they  reported  a  similar 
nonlinearity  (softening)  in  the  pressure  data  as  temperature  is  increased,  and  at  high 
temperature  in  the  range  of  260  °C  and  extending  to  405  °C.  They  reported  that 
their  analysis  of  the  pressure  data  for  SbF  yielded  AHyapor  =  -64  kJ/mol  near  the 
boiling  point  of  401.6  °C. 

Gadzhi  and  coworkers6  measured  the  vapor  pressure  at  a  similar  but  hotter  range 
than  that  investigated  by  Bruner  and  Corbett.  Because  their  data  include  the  highest 
temperatures  ever  recorded,  there  is  likely  some  chemical  decomposition  of  SM3 
(cracking)  into  distinct  species  such  as  free  I  and  Sbl(3-X)  with  a  concomitant  supra 
nonlinear  increase  in  the  vapor  pressure  at  the  highest  temperatures. 

Ferro  and  colleagues7  studied  vapor  pressure  as  a  step  toward  understanding  vapor 
transport  of  chemical  species.  In  contrast  to  the  earlier  work  of  Bruner  and  Corbett, 
Ferro  and  colleagues’  work  focused  on  colder  SM3  solid  near  the  melting  point,  and 
thus  a  much  lower-temperature  range:  starting  at  109  °C  and  extending  to  277  °C. 
Using  2  different  techniques  to  determine  the  equilibrium  vapor  pressure,  they 
reported  very  linear  behavior  and  AHvapor  =  -109.3  kJ/mol . 

The  work  of  Gregory8  did  not  focus  on  the  determination  of  AHyapor-  but  instead 
attempted  to  reconcile  some  of  the  apparent  inconsistencies  in  vapor  pressure  across 
the  entire  temperature  range  of  past  work.  Gregory  used  absorbance  by  the  vapor 
phase  as  an  index  of  vapor  pressure.  One  finding  Gregory  validated  was  that  the 
thermodynamic  behavior  of  SM3  follows  the  trend  of  progressively  reducing 
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AH  vapor  as  temperature  increases.  The  measurements  of  Gregory  were  generally 
consistent  with  past  work,  and  that  suggests  the  nonlinear,  softening  AHyapor  with 
temperature  of  previous  reports  is  generally  consistent. 

5.  Significance 

This  work  expands  the  breadth  of  vapor  pressure  measurement  from  the  previous  6 
orders  of  magnitude  to  16  orders  of  magnitude,  and  represents  a  more 
comprehensive  perspective  of  the  equilibrium  vapor  pressure  of  this  interesting  and 
technologically  important  material.  The  present  work  is  consistent  with  previous 
determinations  of  AHvapor  on  the  order  of  approximately  100  kJ/mol  (Fig.  3).  The 
broadly  arcing  curve  suggests  that  the  AHyapor  value  one  obtains  is,  in  part, 
determined  by  the  temperature  where  the  data  were  measured:  AHvapor  will  be 
higher  at  lower  temperature  and  lower  at  higher  temperature.  Ignoring  entropy 
contributions  to  sublimation  allows  the  contributions  to  the  sublimation  to  be  one 
lumped  term  called  the  effective  enthalpy  of  sublimation.  In  that  case,  the  effective 
AHyapor  of  SbFs  has  a  strong,  nonlinear  temperature  dependence. 
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Fig.  3  Effective  AHvapor  of  Sbl3  as  a  function  of  temperature 
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6.  MBE  Doping  Strategy 


To  be  effective  as  a  dopant,  I  must  be  incorporated  within  the  (Bi,Sb)2(Te,Se)3 
crystal  lattice  during  growth.  If  one  assumes  unity  incorporation  (i.e.,  one  I 
becomes  incorporated  from  one  SbB  molecule  that  hits  the  surface),  then  a  simple 
ratio  of  pressures  could  yield  an  estimate  of  the  obtained  doping  concentration. 
During  (Bi,Sb)2(Te,Se)3  MBE  growth,  the  growth  rate  is  governed  by  the  pressure 
of  Bi  metal.  As  a  consequence  of  that,  the  ratio  of  SbB  pressure  to  Bi  pressure  can 
therefore  be  an  approximate  index  for  the  expected  doping  concentration. 

If,  say,  the  overall  concentration  of  desired  doping  for  I  atoms  is  approximately 
5xl019  I  atoms  per  cm3  and  assuming  each  I  donates  one  extra  electron,  then  that 
would  yield  an  n-type  carrier  density  of  5xl019/cm3.  If  the  (Bi,Sb)2(Te,Se)3  matrix 
concentration  is  roughly  5xl023  atoms  per  cm3,  then  the  preferred  ratio  of  Bi 
pressure  to  SbB  pressure  would  be  10,000:1.  For  a  common  Bi  pressure  during 
(Bi,Sb)2(Te,Se)3  MBE  growth,  a  Bi  pressure  of  5x10 7  Torr  is  not  unreasonable.  In 
the  case  where  Bi  pressure  is  5x10  7  Torr,  one  would  target  an  S b I3  pressure  of 
5x10"' 1  Torr;  refer  to  Fig.  2.  One  would,  therefore,  select  an  SbB  temperature  of 
roughly  -10  °C  to  achieve  an  SbB  pressure  of  5x10  11  Torr,  resulting  in  a  doping 
concentration  of  5xl019/cm3. 

7.  Conclusions 


The  equilibrium  vapor  pressure  of  SbB  has  been  measured  over  a  temperature  range 
from  253  to  320  K  (-20  to  47  °C)  using  a  classical  Knudsen  cell  method.  An 
analysis  of  the  slope  of  that  data  as  a  function  of  temperature  is  one  measure  of  the 
enthalpy  of  sublimation  (AH vapor),  and  a  value  of  AHvapor  =  -128.48  kJ/mol  has 
been  obtained  for  SbB. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


Bi 

I 

MBE 

OMVPE 

P 

p-BN 

Sb 

Sbl3 

Se 

T 

Te 

AH  vapor 

AV 


bismuth 

iodine 

molecular  beam  epitaxy 
organometallic  vapor  phase  epitaxy 
pressure 

pyrolytic  boron  nitride 

antimony 

antimony  triiodide 

selenium 

temperature 

tellurium 

enthalpy  of  vaporization 

net  change  in  molar- specific  volume 
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